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Abstract. 

A unified treatment of all known types of singularities for flat, isotropic and 
homogeneous spacetimes in the framework of loop quantum cosmology (LQC) is 
presented. These include bangs, crunches and all future singularities. Using 
effective spacctimc description we perform a model independent general analysis of 
the properties of curvature, behavior of geodesies and strength of singularities. For 
illustration purposes a phenomenological model based analysis is also performed. We 
show that all values of the scale factor at which a strong singularity may occur are 
excluded from the effective loop quantum spacetime. Further, if the evolution leads 
to either a vanishing or divergent scale factor then the loop quantum universe is 
asymptotically deSitter in that regime. We also show that there exist a class of sudden 
extremal events, which includes a recently discussed possibility, for which the curvature 
or its derivatives will always diverge. Such events however turn out to be harmless 
weak curvature singularities beyond which geodesies can be extended. Our results 
point towards a generic resolution of physical singularities in LQC. 
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1. Introduction 

Singularities are common in General Relativity (GR). These are the boundaries of 
spacetime which can be reached by an observer in a finite proper time where the 
spacetime curvature and tidal forces become infinite. Singularity theorems of Penrose, 
Hawking and Geroch show that the primary characteristic of a physical singularity is 
the inextendibility of the geodesies beyond it [TJ [2] . However behavior of geodesies is 
insufficient to capture the detailed features of singularities and distinguish physical from 
unphysical ones. Hence singularities are also classified in terms of strong and weak types 
0, HI E]. A singularity is strong if the tidal forces cause complete destruction of objects 
irrespective of their physical characteristics, whereas a singularity is considered weak 
if tidal forces are not strong enough to forbid passage of objects or detectors. Due to 
this reason only strong singularities are generally considered as physical. An example 
of strong singularity is the big bang singularity in cosmological models and an example 
of a weak singularity is the shell crossing singularity in gravitational collapse scenarios 
where even though curvature invariants diverge, 'strong detectors' can pass the extremal 
event [6]. 

Traditional cosmological singularities, such as big bang and big crunch, come with 
certain signatures: events where the scale factor goes to zero, geodesies are incomplete 
and objects are crushed to zero volume by infinite gravitational curvature. However, 
recently various new singularities have been discovered in GR [TJ isjji These singularities 
which are typically investigated in the future evolution of the universe (hence popularly 
known as future singularities), do not occur at a vanishing scale factor. The latter 
either goes to infinity in finite proper time, along with a similar behavior of energy 
density p and pressure P (the Big Rip) or the singularity is sudden i.e. at a finite value 
of time and scale factor, curvature or one of its higher derivative blow up [8j [UJ ITT] . 
Features of cosmological singularities can be classified using the triplet of variables 
(a, p, P) and have been understood in detail [121 [13]. For a universe with a Robertson- 
Walker metric and matter satisfying a non-dissipative cosmological equation of state: 
P = P(p), cosmological singularities apart from the big bang and big crunch can be 
completely classified in four types [H]: Type I as big rip, type II as the sudden one 
where energy density is finite but the pressure diverges at the extremal event, type III 
where both energy density and pressure diverge at a finite value of scale factor and 
rest of the remaining as type IV where curvature components are finite but their higher 
derivatives blow up. 

An open question is the way quantum gravitational effects change the picture 
near above extremal events. This can be asked in different ways: Does quantum 
gravity resolve all spacelike singularities? Do quantum gravity effects always bound 
the spacetime curvature? Are geodesies complete (if their notion exists in a quantum 
spacetime)? What does the singularity resolution or the lack of it tell us about the 
underlying theory? Since we do not yet have a complete theory of quantum gravity 

f Similar solutions have also been found in braneworld models [9]. 
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these questions cannot yet be answered in full generality. However, they can be posed 
for cosmological singularities in a simplified setting such as a homogeneous universe 
where mini-superspace quantizations are available. 

Loop quantum cosmology (LQC) [15] is one of the settings where such questions 
can be answered. It is a non-perturbative canonical quantization of homogeneous and 
isotropic spacetimes based on loop quantum gravity (LQG) [16]. The classical phase 
space variables are the Ashtekar connection and the conjugate triad. The elementary 
variables used for quantization are the holonomies of the connection and fluxes of the 
triad. The quantization follows the Dirac's program for constrained systems where 
a physical Hilbert space is obtained by solving all of the constraints at the quantum 
level and predictions are extracted via Dirac observables. In LQC, due to underlying 
symmetries of the homogeneity and isotropy the only non-trivial constraint is the 
Hamiltonian constraint. This is expressed in terms of holonomies and fluxes and then 
is quantized. The quantization has been successfully and rigorously completed in various 
interesting cases which include spatially flat FRW models with at least one massless 
scalar [171 [HI [19], with and without cosmological constant [201 [21], closed [221 123] 
and open universes [21] as well as the inflationary spacetimes [25J. All these models 
classically exhibit big bang (and in some cases, also big crunch) singularity which 
is generically resolved in LQC. As an example, backward evolution of states which 
correspond to an expanding macroscopic universe at late times results in a quantum 
bounce to a contracting universe when energy density becomes equal to a critical value 
p crit ?z 0.41ppi. Various results have been shown to be quite robust using an exactly 
solvable model in LQC (sLQC) with massless scalar as the matter content |26j. In 
particular, one can show that the bounce is a generic property of states in the physical 
Hilbert space and that p crit is the supremum of the spectrum of the energy density 
operator on the physical Hilbert spacqjj The universe in the branch preceding ours has 
been shown to retain semi-classicality [28] and the quantization turns out to be unique 
if one demands consistency and physical viability of the quantization scheme |29j . 

An interesting feature of loop quantum cosmology is the availability of an effective 
spacetime description [301 [31]. Numerical simulations of the exact LQC equations for 
universes which grow macroscopic at late times confirm that the effective modified 
Friedman dynamics captures the underlying quantum dynamics and the bounce to an 
amazing accuracy [HI [191 [231 [201 Ell 121] • These results were obtained using massless 
scalar field with and without cosmological constant and also the inflationary potential. 
Assuming that key features of the effective dynamics can be trusted more generally, 
various interesting results have been obtained. These include the resolution of big 
crunch singularities for negative potentials in Cyclic models [32] as well as in string 
inspired pre big bang scenarios [33] and relation with Palatini theories [34|. 

In the dark energy scenarios, phantom field models have been analyzed at an 
effective level and are shown to be generically free of big rip or the type I singularity 

| The boundedness of the density operator also holds for a wide class of lapse function [27] . 
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[3l)| EE] • An interesting example which is studied in this context is the case of a phantom 
model with an unbounded negative (positive) potential for a canonical (phantom) field 
[37] . Classical dynamics for this model predicts a big rip singularity in future. Using 
effective dynamics of LQC, authors of Ref. [37J find that as the future singularity is 
approached, though the energy density is bounded, the pressure and the rate of change of 
the Hubble rate blow up in LQC and they conclude that the "singularity" is unavoidable. 
This case is interesting because it is an explicit example where the curvature invariant 
is not bounded by the quantum geometric effects. For the particular model considered 
in Ref. [37], quantum geometric effects convert the classical type I singularity in to 
a sudden type II singularity. As we will show, the latter "singularity" is weak and 
unphysical. 

Although the fate of singularities has been studied for specific models in LQC, 
a general treatment has so far been unavailable. Further, details about the nature 
and strength of the various singularities and the properties of geodesies in the effective 
spacetime were so far not investigated. Given that an effective spacetime description 
is available in LQC, all these issues can be addressed and analyzed in detail. This is 
not only important to distinguish physical singularities from unphysical ones but also 
to understand the way loop quantization affects the fate of spacetime at these singular 
events. Moreover the examples which have so far been studied do not fully exhaust all 
the possibilities which include type III and type IV singularities. The aim of this work 
is to understand all these issues for all cosmological singularities in flat (k=0) model of 
LQC with matter satisfying a non-dissipative cosmological equation of state P = P(p). 
Our work will assume (i) that the effective spacetime description is valid for all matter 
staisfying above equation of state and (ii) the effective value of geometrical and physical 
entities such as geodesies and curvature invariants coincide with those derived from the 
effective spacetime metric jjj] 

We organize this paper as follows. In the next section we briefly revisit loop 
quantization and the effective Friedman dynamics in LQC (for details we refer the reader 
to Ref. [29J ) . In Sec. Ill we review all the singularities in k = FRW model and derive 
the geodesic equations for the Robertson- Walker metric and state the Clarke-Krolak 
conditions [38] which are necessary as well sufficient for a singularity to be considered 
strong a la Tipler [I] and Krolak [5]. In Sec. IV we consider a model which is sufficiently 
general to illustrate the singularity resolution in LQC. This model was introduced in the 
context of dark energy scenarios [2] and exhibits all possible cosmological singularities 
including both strong and weak. We then show that in this model all strong singularities 
are resolved. The weak singularities however remain unaffected. In Sec. V we provide 
a model independent analysis using effective dynamics of LQC. Here we first show 
that in the effective spacetime of LQC, if the evolution is such that the scale factor 

|| In the models with non- vanishing intrinsic curvature, the following analysis will also require 
implementation of inverse scale factor effects which modify the matter conservation equation and may 
become dominant when scale factor approaches zero. In the flat model, unless one restricts to a compact 
topology these effects are argued to be unphysical [19] and are not considered here. 
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either vanishes or becomes infinite then the universe always approaches a deSitter state. 
Cosmological observers in an effective loop quantum spacetime take infinite proper time 
to reach above values of the scale factor. As in GR, these are nonsingular. In most cases 
loop quantum evolution does not lead to an asymptotic deSitter regime and these are 
carefully analyzed using Lipshitz conditions for dynamical as well geodesic equations. 
Using the fact that energy density and hence Hubble rate are always bounded above in 
LQC, it is shown that dynamical and geodesic equations never break down. For above 
cases geodesies can be extended to arbitrary values of the affine parameter. We show 
that all strong singularities are generically resolved in flat and isotropic LQC. The only 
possible "singularities" are weak curvature type. These are harmless as geodesies can 
be extended beyond them. We summarize the results with a discussion in Sec VI. 



2. Preliminaries 



We will consider the dynamical features of k = homogeneous and isotropic universe in 
LQC. The 3+1 spacetime is described by the manifold Exl, where S is the non-compact 
spatial manifold, and the Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t) (dr 2 + r 2 (d6 2 + sin 2 Odcf) 2 )) (1) 

where a(t) is the scale factor of the universe. Here we have chosen lapse iV to be equal 
to unity so that t is the proper time. 

Effective description for the loop quantization of above spacetime can be obtained 
using geometric formulation of the quantum theory. Using coherent state techniques it 
is possible to derive an effective Hamiltonian (up to controlled higher order corrections) 
for various matter sources [30]. It turns out that for states which correspond to a 
macroscopic universe, such as ours, at late times the following effective Hamiltonian 
captures the underlying loop quantum dynamics: 

-H - 3 Sin2 ^ y I T-f (2) 

Here (3 and volume V = a 3 are conjugate variables satisfying 

{(3,V}=4nG 7 (3) 

with 7 rs 0.2375 as the Barbero-Immirzi parameter. In the classical theory the phase 
space variable /3 — ja/a. The parameter A captures the discreteness of the underlying 
quantum geometry and its value is determined by the minimum eigenvalue of the area 
operator in loop quantum gravity (LQG) [39|, 

A = 2(v / 3vr 7 ) 1 / 2 £ P1 . (4) 

We consider matter to be minimally coupled and homogeneous. In particular it satisfies 
a cosmological equation of state P = P{p) where P is its pressure and p is its energy 
density. 

It is to be noted that though there may exist further state dependent quantum 
corrections to the effective Hamiltonian, the numeric simulations which have so far 
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been performed show that they turn out to be negligible for states representing realistic 
universes we are interested in (see for example Refs. [23J, [20l EH [25] and also refs. [4"0ll4"T] 
for anisotropic models). Guided by these results our analysis will assume the existence 
of above effective Hamiltonian for general matter. 

Given eq.(T5]), it is straight forward to find the modified Friedman dynamics. The 
vanishing of the Hamiltonian constraint Ti, e s ~ leads to 
sin 2 (A/3) 87rG 7 2 

where p = TC mat t/V is the energy density. Then from the Hamilton's equation 
V = {V,H cS } = -AttG^-^Hcs 

3 sm(\3) , , _ , 

= V— - cos(A/3) V (6) 

7 A 

we can obtain the modified Hubble rate 

where p cr i t is given by 

Pcrit = 3/(87rG 7 2 A 2 ) « 0.41p P1 . (8) 

A similar calculation for the second Hamilton's equation: (3 = {P,H e s} results in 
the modified Raychaudhuri equation 

* = p (l - 4JL) - 47 rG P (l - 2-0. (9) 



Pcrit / \ Pcrit , 

Since quantum geometry does not affect the matter part, the Hamilton's equation for 
matter field yield the conservation law 

p + 3H(p + P) = (10) 

where pressure P = —dH matt /dV. It is straight forward to check that Eqs. (j7|), © and 
(flO]) form a closed set. 

The modified Friedman and Raychaudhuri equations are sufficient to determine the 
non-trivial components of the Ricci (and the Einstein tensor) on the FRW background 
with modified expansion rate of the scale factor. The Ricci curvature invariant turns 
out to be 

' r2 , a \ _ o_^„ fi o„.. , o P 



R = 6 [H* + - = 87rGp 1 - 3w + 2^— (1 + 3w) . (11) 

V <V V Pcrit / 

where w is the equation of state of the matter component w = P/p. 

Classical equations for the FRW spacetime can be obtained from Eqs. (1Tf9"j) and (TTTT) 
in the limit A (that is Gh — > 0): 

H 2 = —p, - = - p + 3P 12 
3 a 3 

and 

R = 8nG(p - 3P) . (13) 
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The effective dynamical equations immediately lead to an upper bound on the 
energy density and hence the Hubble rate in LQC As we will show in Sec. IV and Sec. 
V, these features play an important role in absence of physical singularities in a loop 
quantum universe. 

3. Cosmological Singularities: Nature and Strength 

Singularities for the homogeneous and isotropic spacetime can be classified using the 
behavior of scale factor, energy density and pressure (or equivalently in terms of 
spacetime curvature). All of the known (and plausible) singularities for matter satisfying 
non-dissipative equation of state P = P(p) fall in one of the categories below [T2l [TBI [hi] : 
Big Bang and Big Crunch: These are accompanied by vanishing of the scale factor 
at a finite proper time and the divergence in energy density and curvature invariants. 
Null energy condition (NEC), (p + P) > 0, is always satisfied in these events. 

Big Rip or Type I singularity: For these singularities NEC and hence all other en- 
ergy conditions are violated. The scale factor diverges in proper finite time, a(t) — ► oo. 
This is accompanied with a divergence of energy density, pressure and curvature invari- 
ants. 

Sudden or Type II singularity: This extremal event occurs at a finite value of the 
scale factor a — > a e . It is characterized by a finite value of the energy density but an 
associated divergence of pressure. Due to the latter, R diverges. 

Type III singularity: As type II singularity, this singularity also occurs at a finite 
value of the scale factor. However, both the energy density and pressure diverge, causing 
a blow up of curvature invariants. 

Type IV singularity: None of the energy density or pressure blow up in this case 
which occurs at a finite value of the scale factor. Curvature invariants are finite, however 
curvature derivatives blow up. 

Singularities associated with a divergence in spacetime curvature fall in either big 
bang/crunch or type I-III class. Remaining singularities are curvature derivative kind 
which form Type IV class. Unlike big bang/crunch and type I singularities, analysis of 
energy conditions for type II, III and IV is subtle and answers can be model dependent 
[12] . However, it is always true that Type II singularities are accompanied by violation 
of dominant energy condition (DEC): (p ± P) > 0. 

To understand the behavior of geodesies let us consider the geodesic equations for 
the flat (k = 0) Robert son- Walker metric: 



where prime denotes a derivative with respect to the affine parameter (r). Using 




(14) 
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Cartesian coordinates, u x ,u y and u z turn out to be constant. Therefore it is sufficient 
to analyze geodesic equation for time coordinate which is 

** = ° + W) (15) 

where x is a constant and e = 1 for massive particles and e = for null geodesies. Since 
below we also consider radial geodesies, it is also useful to obtain equation for radial 
geodesies 

a 2 (t) r' = X (16) 
which using (IT51) implies 

f = -aar' 2 = -H (t' 2 - e) . (17) 

Understanding properties of above geodesic equations is important to prove whether 
the spacetime is geodesically complete. For that it is necessary to show that geodesic 
equations admit a unique extendible solution. From the geodesic equations we see that 
these break down when scale factor becomes zero and/or the Hubble rate blows up. 
Hence these break down at big bang/crunch, type I and type III singularities. Since 
type II and type IV singularities occur at a finite value of scale factor and Hubble rate 
remains finite at these events, therefore geodesic equations do not break down. We will 
later show that for these cases a unique extendible solution to the past and future of 
type II and type IV singularities can be found (see also Ref. [12] )■ Thus, type II and 
type IV singular events fail to be singularities a la theorems of Penrose, Hawking and 
Geroch. From the criteria of geodesic inextendibility only big bang/crunch, type I and 
type III singular events turn out to be singularities. 

Apart from analysis of geodesies it is also useful to consider the strength of the 
singularities. A detailed discussion of the strength of various cosmological singularities 
is available in Ref. [13] . We here only note the conditions necessary for our analysis. 
These originate from the work of Clarke and Krolak to differentiate types of singularities 
which involve analysis of integrals of curvature components for both null and particle 
geodesies [38]. For the FRW metric, a singularity occurring at the value of the affme 
parameter r = r e is a strong curvature type a la Tipler iff the following integral over 
the spatial components of Ricci tensor is unbounded 



air 1 / dr"R ab u a u b , (18) 
'o Jo 

as r — > r e . Else the singularity is weak. A less restrictive condition is by Krolak who 
classifies the singularity to be strong iff 

drR ab u a u b (19) 



o 



is infinite as r — > r e . It is thus possible that a singularity may be strong a la Krolak 
but weak a la Tipler. A strong singularity from above conditions is the one in which an 
in-falling observer or detector is completely annihilated by the tidal forces. For a weak 
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curvature singularity, tidal forces are not strong enough to cause such a destruction. 
Sufficiently strong detectors survive such events. 

Note that the integrand in both of above integrals is proportional to combination 
of square of Hubble rate and d/a. As an example, for null geodesies 

^^W)( H *-%) (20) 

where we have used (IT51) and ffTEj) . Since conditions (fT8|) and (TTTJ|) involve at least one 
integral over affine parameter, it turns out that the integrals are finite if only a or higher 
derivatives diverge and the scale factor neither vanishes nor diverges. This happens in 
the case of type II and type IV singularities which are thus weak. Big bang, big crunch 
and big rip are strong curvature singularities according to both Tipler and Krolak. 
Type III singularities are strong according to Krolak's condition but weak by Tipler's 
condition [13]. (Similar conclusions are reached by the analysis of particle geodesies). 
It is to be noted that strong curvature singularities in FRW universe are also the ones 
beyond which geodesies can not be extended. On the other hand weak singularities are 
the ones beyond which geodesies can be extended and thus are harmless events. 



4. Singularity resolution in LQC: Illustration via a Model 

We now illustrate loop quantum dynamics and results of the previous section using 
a general model which exhibits all the cosmological singularities of interest |14j . This 
model is of interest because it describes a general dark energy scenario in a FRW universe 
including quintessence and phantom dark energy models. Classically the model predicts 
various singularities for different ranges of parameters and it proves useful to understand 
the detailed properties and the fate of the universe in such scenarios. The model is based 
on the ansatz 

P = -p-f(p) (21) 

with 

'M = ■ < 22 > 

Here A, B and a are parameters of the model. Their values determine the nature of 
singularities. Note that when f(p) = 0, the model reduces to the standard cosmological 
constant scenario. Above anasatz therefore proves very useful to study departures of 
the equation of state from the cosmological constant setting. 

The dependence of energy density on scale factor can be found by integrating ffTUl) 

a = a Q exp -— f f (23) 



which yields 



6 AB(l-a) 



l/2\ 
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To investigate resolution of various singularities we will use eqs.flT]), (ITT]) along 
with R = 6(H + AHH) and 

P = ±v / 247rGp(l-p/p crit )(p+P) 

In the following we will only consider the future evolution and compare the results of the 
classical and loop quantum evolution. (For details of the classical dynamics for various 
choices of parameters of this model we refer the reader to Ref. p3].). 



{za — i ) /\a pr~ 
Art*- 1 + B 



\l — a.)A np 



4-1. Type I Singularities 

If the value of a is chosen between 3/4 < a < 1 and A is positive, then the model gives a 
big rip (type I) singularity in GR. The scale factor, energy density and pressure diverge 
at a finite time and the DEC is violated (for all times). There is no big bang in the 
classical theory (since DEC is violated). The model is devoid of an initial singularity. 

In LQC, the big rip singularity is avoided. The energy density initially grows as 
in the classical theory, however when it becomes comparable to p C rit ; departures from 
classical trajectories become significant. Eventually, p becomes equal to p crit and the 
Hubble rate vanishes with a taking negative value. The universe instead of ripping apart 
in finite time, recollapses and the evolution continues. The Ricci scalar, its derivatives 
and higher curvature invariants are bounded in the entire evolution. 

In Fig. [T] we compare the evolution of Hubble rate and Ricci curvature scalar in 
the classical and the effective dynamics of LQC. The classical Hubble rate diverges as 
a — > oo, whereas the Hubble rate in LQC is bounded. As can be seen, unlike in GR, R 
is bounded in the loop cosmological evolution. 



4-2. Type II Singularities 

A necessary condition for type II singularities to occur is A/B < 0. For these cases 
there is also a big bang (crunch) singularity as a — > 0. The sudden singularity occurs 



Are loop quantum cosmos never singular? 



11 



at a finite value of the scale factor, a —* a Q . As the singularity is approached the energy 
density goes to zero however pressure and hence the Ricci curvature diverge. Since the 
Hubble rate is bounded, the geodesies are extendible and the singularity is only a weak 
curvature singularity. 

In loop quantum evolution the initial singularity which is a strong curvature type is 
resolved. Since p <C p crit when the sudden singularity occurs, the latter is not resolved. 
In fact near this extremal event the dynamics mimics the one obtained from GR and the 
properties of geodesies do not change qualitatively near the sudden singularity. As in 
the classical theory, a/ a and R diverge near the type II event. Since the initial big bang 
singularity is also resolved the effective spacetime in LQC is geodesically complete. 

Dynamical trajectories obtained from numerical integration of classical and loop 
quantum equations are compared in Fig. [2j The classical Hubble rate (dashed) diverges 
at the big bang and approaches zero near a Q . The loop quantum Hubble rate (solid 
curve) is bounded throughout the evolution and agrees with the classical values near 
a — ► a Q . It vanishes at the bounce point in the early universe and approaches the 
classical value at late times. The behavior of Ricci scalar shows that it is bounded in 
LQC in the early epoch and attains a positive value at the quantum bounce. However, 
it diverges when the sudden singularity is approached. The quantum geometric effects 
are unable to bind the curvature in this case. 



3.0 ; S n * 

2-5: ""-^ 100- 




2000 4000 6000 8000 10000 

a -200 L 



Figure 2. Plot of Hubble rate and Ricci scalar for type II singularity. The classical 
Hubble rate diverges at the big bang but is finite (goes to zero) at a = a a . Parameters 
are A = -0.1, B = 1 and a = 1/4.1. 



Since pressure grows unboundedly as the type II singularity is approached, there 
is a huge violation of DEC near a . For a < 0, using (T2T1) we obtain: w — ► — oo for 
positive B when a — > a Q . It can also be shown that these extremal events do not satisfy 
Tipler and Krolak's conditions for strong singularities. For null geodesies we obtain 

RiM = 8nG(p + P)^r(l- 2-0 , (26) 

a 2 V Pcrit/ 

using which we can compute the Tipler [I] and Krolak [5] integrals. Both of these are 
finite in LQC. As an example we show the integrand of ( fT9l) (after a change of variables 
to the scale factor) in Fig. [3j Value obtained from numerical integration for the chosen 
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values of parameters turns out to be approximately 5.6 x 10 4 . Similar results follow 
for the particle geodesies. 



H i.o 



60 



R 




Figure 4. The classical (dashed) Hubble rate and Ricci scalar is compared with the 
one in LQC for type III extremal events. Classically Hubble rate blows up as the 
singularity is approached at late times. Both the Hubble rate and Ricci are bounded 
in LQC. Parameters are A = 100, B = 1 and a = 2. 



4-3. Type III Singularities 

When a > 1, the model predicts a type III singularity at a finite time and scale factor 
a — > a a . Both the energy density and pressure diverge in GR. The unbounded Hubble 
rate causes geodesies to be incomplete and results in a strong curvature singularity a la 
Krolak. Since the singularity occurs at a finite volume, it is weak according to Tipler's 
criteria. In this model the initial singularity is absent at the classical level itself. 

In loop quantum dynamics there is no type III singularity. When p is small, the 
classical and effective evolutions are similar. However they become qualitatively different 
as p increases. The loop quantum universe recollapses at p — p cvl t and the type III 
singularity is avoided. The Hubble rate, a /a and Ricci scalar are bounded and finite. 
Thus geodesies are complete and integrals ([18]) and f|T9|) are finite. 
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A comparison of the classical and effective LQC evolution is presented in Fig. HI 
As can be seen in the plot the Hubble rate is bounded in LQC, whereas it grows un- 
bounded classically in the future evolution. Similarly, the Ricci scalar which diverges in 
the classical theory clS Qiq IS approached is bounded in LQC. 



4-4- Type IV Singularities 

If the value of a is between < a < 1/2, then though the energy density 

and pressure remain finite, a higher derivative of curvature diverges. This singularity is 
a derivative curvature singularity and none of the curvature invariants diverge. As in 
type II case, since Hubble rate is finite at a = a a , geodesic equations are well behaved. 
The singularity is a weak curvature type under the classifications of both Tipler and 
Krolak. 

3.0 ; \ 200 - 



2.0 1 




Figure 5. Comparison of Hubble rates and R is made for the classical theory and LQC. 
The classical divergence is controlled by the quantum geometry leading to a bounce in 
the early universe. At late times, the classical curve is a good approximation to LQC 
and both vanish at a = a a . Parameters are A = —0.1, B = — 1 and a = 1/4.1 



Quantum geometric effects have little influence on this harmless extremal event 
beyond which geodesies can be extended even in the classical theory. However they do 
resolve the big bang singularity for this model which accompanies a classical divergence 
of Hubble rate and Ricci curvature and thus lead to a geodesically complete spacetime. 
As is evident from the Fig. [5j Hubble rate is finite through out the loop quantum 
evolution. 

For type IV singularities, the value of a determines the order of derivative which 
blows up at a = a D [T4J. For a = 1/4.1, the divergence occurs in H and hence R which 
is depicted in the second plot of Fig [5] 

5. Singularity resolution in LQC: General analysis and some observations 

In the previous section we saw that the effective dynamics of LQC successfully resolved 
all types of strong singularities in FRW cosmology for quite a general equation of state 
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in a model proposed in Ref. [H]. An important question is whether these results hold 
in general for a cosmological equation of state of the form P = P(p). To understand 
this let us first note that modified Friedman dynamics in LQC (J7|) leads to a universal 
upper bound p CTit , i.e. 

< p < Pcrit (27) 

and the Hubble rate is bounded with the maximum allowed value at p = p cr \ t /2: 

1/2 



U = \ v/316^ 7 3 ) (28) 
(where we have used the expression for p crit ([8])). This implies a bound on the trace 
of the extrinsic curvature K, since for the flat Robert son- Walker metric it is related to 
Hubble rate as K = 3H. Note that the upper bound arises purely because of quantum 
gravity. If Gh —>■ we find that |-£/"| m ax — * oo. 

Since energy density is bounded above by p crit , all values of the scale factor for 
which p > p crit are excluded from the effective spacetime of LQC. This immediately 
implies that singularities associated with divergence of energy density and Hubble rate 
are absent in LQC. These include big bang/crunch, type I and type III singularities i.e. 
all strong singularities in FRW cosmology. 

Before we analyze the nature of geodesies in LQC, let us note an interesting property 
of the effective equations related to two values of the scale factor: a(t) = and a(t) = oo. 
In classical cosmology depending on the equation of state of matter, a big bang/ crunch 
singularity may occur at a(t) = and a big rip (type I) singularity may occur at 
a(t) = oo. However, in effective loop quantum spacetime this is not the case. In LQC, 
if evolution leads to above values of scale factor then the universe always approaches a 
deSitter state, i.e. equation of state becomes that of the positive cosmological constant 
(w = —1). From the modified Friedman dynamics we find that unlike classical theory, 
energy density, Hubble rate and curvature invariants are always finite at a(t) — and 
a(t) = oo in LQC. It can also be shown that cosmological observers take infinite proper 
time to reach above values of scale factor in the effective loop quantum spacetime. 

We prove this by contradiction. Let us first note that the following identity obtained 
from the conservation law ffTOl) holds for all values of the scale factor: 



\n(-t) =-3 f\l + w (a))^ . (29) 
\PoJ J ao a 

We assume that w(a) be a smooth real function of the scale factor. Let us now assume 
that it is possible to approach a(t) = without w — > — 1 and satisfy above equation. 
Since left hand side of above equation is finite, it requires the integrand to be such that 
the right hand side is also finite. This implies that the product of above integrand with 
scale factor should go to zero as a — > 0. Which means that (l + w(a)) must vanish in the 
above limit for Eq. (l29|) to be satisfied, i.e. w — >• — 1 as a(t) — » 0. If w ^ —1 as a(t) — > 0, 
the right hand side blows up and the equation is not satisfied. Hence our assumption is 
incorrect. 



Are loop quantum cosmos never singular? 



15 



Similarly we can prove that in LQC it is not possible to approach a = oo without 
w(a) — > — 1. To prove it by contradiction, let us assume that we can approach a(t) = oo 
without w — > — 1 and satisfy eq.( l29i) . Since the R.H.S of ( 1291) is finite, the L.H.S is 
finite only if the product of the scale factor and the integrand goes to zero as a — > oo. 
That is, (1 + w(a)) — > as a(t) — > oo. For equation of state not approaching —1 as 
a — > oo, above equation can not be satisfied and therefore our assumption turns out to 
be incorrect. Hence we can state, 

Remark 1 In flat isotropic LQC if the evolution leads to either a vanishing or a 
divergent value of the scale factor then the universe is asymptotically deSitter in that 
regime. 

This feature of LQC is very interesting. To understand it futher let us consider loop 
quantum dynamics with a positive cosmological constant (A). Here we first note that 
for a matter as pure cosmological constant the modified Friedman and Raychaudhuri 
equations in LQC are equivalent to the classical ones with a "renormalized" cosmological 
constant. To see this let us consider the modified Friedman and Raychaudhuri equations 
for a deSitter universe sourced with energy density p = A/(8irG), 

9 A' , a A' , s 

H 2 = - and - = - 30 
3 a 3 

where A' is the "renormalized" cosmological constant^ 

A' — a(i — — — ^ ) . (31) 



v 8nGp crit/ 

Solving above field equations we can obtain the solution to the scale factor in the small 
neighbourhood of a(t) = or a(t) = oo. The solution behaves as a(t) ss exp(±-^/A'/3 t) 
as we approach above values of the scale factor. Hence cosmological observers in LQC 
take infinite proper time to reach a(t) = or a(t) = oo. As in classical GR, in these 
cases the spacetime is extendible (and in this sense it is non-singular). 

Above cases where the evolution leads to a vanishing or a divergent scale factor are 
not so common in a loop quantum evolution. In most cases of interest evolution does not 
lead to above values and we now focus on them. To understand the behavior of geodesies 
and their extendibility in these cases, it is useful to first consider the dynamical equations 
(J7|) and (flUj) and analyze Lipshitz conditions for existence of a unique solution. It is 
straightforward to find that except the following critical points, the Lipshitz conditions 
are always satisfied and equations are regular. These points are (i) when energy density 
becomes equal to the critical energy density p crit , (ii) when pressure P diverges with 
a finite value for energy density and (iii) when P diverges at a finite value of energy 
density and pressure. First critical point corresponds to the bounce/recollapse point in 
an LQC universe. From the classification in Sec. Ill we see that the second and the 

+ The "renormalization" of cosmological constant due to quantum geometry effects has an interesting 
feature. If A < %/3/(327r 2 7 2 ^ 

|j then the correction to cosmological constant is very small. However, 
if A is of the order Planck or more precisely A w \/3/ (32n 2 "/ 2 £p l ) then the "renormalized" value of 
cosmological constant is very small compared to Planck scale. 
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third critical points correspond to type II and type IV singularities respectively. It is 
important to note that none of the critical points correspond to a strong singularity. 
Further above critical points are not problematic. In the neighbourhood of these points 
we can use 



to determine the scale factor and obtain a unique solution a(t) in the past and future of 
the critical points. Since Lipshitz conditions are satisfied everywhere else except above 
points it is hence possible to obtain a global solution for the dynamical equations. 

Now we analyze geodesic equations. The upper bound on the Hubble rate ensures 
that these never break down in LQC. Lipshitz conditions are satisfied for null and 
particle geodesies and hence a unique extendible solution exists. As an example, let 
us consider equations for null geodesies for time ( fl5|) and radial coordinates ([IB]) : 
t' = x/ a (t) = f(ti T ) an d t' = xjd^it) = g(r,r). For cases under consideration: 
f,g < oo. The derivative of / and g with respect to radial coordinate is trivially 
zero. Due to boundedness of Hubble rate, the derivatives with respect to time are finite 
and thus Lipshitz conditions are satisfied. A similar analysis can be performed for the 
particle geodesies and we obtain the same result. Therefore for cases under consideration 
geodesies can be extended to arbitrary values of the affine parameter. 

These results should be contrasted with those in the classical theory where geodesies 
can be extended only beyond type II and type IV singularities, i.e. weak singularities 
P3J H2]. Since geodesies can not be extended beyond strong singularities which 
are common in the classical theory, classical spacetimes are in general geodesically 
incomplete. 

Let us now consider the behavior of curvature invariants. Due to underlying 
symmetries of the Robert son- Walker metric, it is sufficient to analyze Ricci curvature 
scalar (R). As we have seen Hubble rate is always bounded in LQC. If a /a is also 
bounded, then R is bounded. From Eq. liTTj) it is clear that divergence in R can only 
arise if pressure and hence the equation of state diverges. This corresponds to the type 
II singularity in LQC (which was also observed in ref. [3?]). Since geodesies can be 
extended beyond them, these are harmless events. Note that for any reasonable form 
of matter, the equation of state is expected to be finite and hence for such reasonable 
forms of matter curvature invariants never diverge in LQC. 

Finally let us turn to the strength of these extremal events. From discussion of 
Tipler ffl8l) and Krolak ffT9]) conditions in Sec. Ill we know that existence of a strong 
curvature singularity requires a divergence in Hubble and/or a vanishing of the scale 
factor. In LQC Hubble rate is always bounded. The scale factor vanishes or diverges 
only when universe approaches a deSitter phase for which the integrals (fl8|) and (fl9l) are 
finite. Thus Tipler and Krolak integrals are always finite in LQC. Hence we conclude, 




(32) 



and 




(33) 
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Remark 2 No strong curvature singularities exist in the effective spacetime of flat 
isotropic LQC. 

It is interesting to note that loop quantum modifications only resolve strong 
singularities in cosmology whereas weak singularities may still occur in the effective 
spacetime. However as we have shown geodesies can be extended beyond the latter 
events and sufficiently strong detectors survive tidal forces. Thus quantum geometry 
is able to distinguish between physically relevant strong singularities from unphysical 
weak ones and resolves only the former. 

6. Conclusions 

One of the key predictions of LQC is the resolution of big bang (big crunch) singularity 
and existence of quantum bounce at Planck scale in various models pT9] . [20l [2TI [23l [251 [26] . 
The underlying quantum dynamics for states which lead to a macroscopic universe 
at late times can be described by modified Friedman and Raychaudhuri equations 
derived from an effective Hamiltonian [30]. The latter turn out to be very successful 
in capturing details of the true quantum evolution, a feature which has been tested for 
various forms of matter ranging from equation of state of massless scalar to that of the 
cosmo logical constant. Modified Friedman dynamics is hence a valuable tool to probe 
the way quantum geometric effects resolve the singularities in the effective spacetime. 
An important question is whether the results of singularity resolution are generic. The 
aim of this work was to understand this issue in the effective spacetime description of 
flat homogeneous isotropic LQC with a general non-dissipative cosmological equation 
of state. 

Assuming that the effective equations are valid for a general matter model we 
analyzed in detail nature and strength of all possible cosmological singularities and 
behavior of geodesies in the k = isotropic and homogeneous model % Effective 
equations were used to analyze a model with a general enough ansatz for equation 
of state which allows study of all possible singularities in Sec. IV. A general analysis for 
the cosmological equation of state P = P(p) using effective equations was performed 
in Sec. V. We found that singularities which involve divergence of energy density (or 
Hubble rate) are resolved. The underlying reason is the existence of upper bound on 
the energy density of the matter in LQC which translates in to an upper bound for 
the trace of extrinsic curvature. We show that it results in a generic resolution of all 
strong curvature singularities and finiteness of spacetime curvature for flat isotropic 
LQC. Points in the classical spacetime where such singularities occur are excluded in 
the effective spacetime of LQC. However, there do exist extreme events such as type 
II singularities where the curvature diverges due to unbounded pressure but a finite 
energy density. We showed that for such cases quantum geometry plays little role and 



jj Some of these singularities may be already restricted in the full quantum theory due to properties of 
matter Hamiltonian. However to keep the analysis general we allowed all possibilities. 
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divergence in curvature will not be controlled. Similarly, type IV events which involve 
a divergence in curvature derivative can occur in LQC Quantum geometry does not 
exclude these points from the effective spacetime. Interestingly these events are not real 
singularities even in the classical theory and geodesies can be extended beyond them 
[T3l l4"2] . Further, tidal forces are unable to destroy sufficiently strong detectors and 
hence these singularities are weak. 

We also found some nice properties of the effective spacetime in this analysis. First 
that if evolution leads to either a vanishing or a divergent values of the scale factor 
then the loop quantum cosmos behaves as a deSitter universe in that regime. As 
in the classical cosmology, the spacetime can be extended in these cases. Secondly, 
quantum geometry is able to distinguish physically relevant strong singularities from 
weak singularities. Thus, LQC resolves only the real singularities and ignores the 
harmless extremal events. In general loop quantum effects may either completely 
eliminate all strong curvature singularities (as demonstrated by the model in Sec. IV) 
or convert them to harmless weak ones (as in the analysis of Ref. [37]). These results 
are the examples that there can exist physically interesting scenarios where divergences 
in curvature may not be regulated and yet there may be no physical singularity. 

Results obtained in this work can be extended in a straightforward way to the 
curved spatial manifolds (k = 1 models) and Bianchi-I anisotropic spacetimes in LQC 
[43] . It will be interesting to extend the present investigations to models which go 
beyond homogeneous spacetimes and also matter models with a more general equation 
of state than P = P(p). Also, incorporation of higher order state dependent quantum 
corrections to the effective Hamiltonian is a useful direction to explore as they will give 
useful insights on the role of states in singularity resolution [31J. If we look at the way 
loop quantization is performed, the deeper reason behind singularity resolution in flat 
cosmological model is tied to the careful quantization of the Hamiltonian constraint 
in LQC [IH] which turns out to be unique in various ways [29]. Thus leading to a 
harmonious convergence of various results. It is straightforward to see that a different 
choice of quantization which is inequivalent to the improved dynamics [19] (or sLQC 
|26j) would not lead to generic resolution of singularity at the effective level. Hence 
investigations carried here have potential lessons also for the full theory in relation to 
restrictions on various quantization ambiguities. 

We will also like to point out that the bound on the Hubble rate in LQC can 
be viewed as the one on the expansion of congruences in the effective spacetime [44]. 
Whether or not this is a generic feature of a quantum theory of geometry is an 
important open question. Its answer should provide insights on much sought quantum 
generalization of the classical Raychaudhuri equation (Q which is crucial to prove a 
non- singularity theorem |45j . It also remains to be seen whether the simplistic way 
of quantum geometry to resolve the singularities, i.e. ensuring geodesic completeness 
and not a bound on spacetime curvature in general and removing only strong curvature 
singularities, is a common feature of loop quantum spacetimes. In our opinion this 
should serve as one of the guiding principles for such quantizations. 
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